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It is known that quantum fidelity, as a measure of the closeness of two quantum states, is oper- 
ationally equivalent to the minimal overlap of the probability distributions of the two states over 
all possible POVMs; the POVM realizing the minimum is optimal. We consider the ability of ho- 
modyne detection to distinguish two single-mode Gaussian states, and investigate to what extent it 
is optimal in this information-theoretic sense. We completely identify the conditions under which 
homodyne detection makes an optimal distinction between two single-mode Gaussian states of the 
same mean, and show that if the Gaussian states are pure, they are always optimally distinguished. 



PACS numbers: 03.67.Hk, 03.65.Ta, 42.50.Dv, 89.70.+C 



I. INTRODUCTION 

In the field of information processing, both classical 
and quantum, it is crucial to be able to distinguish be- 
tween two items of information. Various distance mea- 
sures have been proposed to quantify the separation of 
items of information, and in particular, the trace distance 
and quantum fidelity are widely used in the quantum in- 
formation field [l| . The quantum fidelity was first defined 
by Jozsa 0], based on Uhlmann's transition probability 



I 



: given two quantum states, p\ and p2, it is given by 
F{p\,p2) = trWi/pi pi -\fp\. The fidelity quantifies the 
closeness of the two states, so that their separation can be 
measured, e.g., by the Bures distance D 2 B = 2(1 — F) Q. 

The quantum fidelity can be given a physical or oper- 
ational significance as follows [a]. By performing a gen- 
eral POVM measurement {E m }, where the E m sum to 
the identity operator [fj, one obtains probability distri- 
butions Pi(m) = t l c(p 1 E m ) and Piijm) = tr(p 2 -£'m) f° r 
the states p± and pi, respectively. It is known that the 
overlap of these distributions is always greater than or 
equal to the fidelity — i.e., J2 m V Pi( m ) V P2{fn) > F. In 
particular, Fuchs and Caves have proved that there al- 
ways exists a POVM satisfying the equality 0, ■ The 
fidelity, then, is equal to the minimal overlap of the de- 
fined probability dis tributi o ns over all possible POVMs, 
F = min {Em} J2 m \fp~Jjn) y/p^m). 

In the proof given by Fuchs and Caves [B| , the optimal 
POVM — that yielding the quantum fidelity — depends on 
the pair of states being compared In experiments, 
however, one would prefer to take a predetermined class 
of measurements, rather than set up the measurement 
case by case. Then, the pairs of quantum states that can 
be optimally distinguished by the chosen measurements 
will be restricted in general, and the condition for such 
pairs depends on the kind of measurement under consid- 
eration. In this work, we consider homodyne detection, 
a continuous variable (CV) measurement, as a means of 
distinguishing two single-mode Gaussian states in the 
aforementioned information-theoretic sense. There has 
been significant progress in quantum information pro- 
cessing using continuous variable systems [jj, and much 



attention has been paid to their Gaussian states which 
are rather easily accessible in experiments and tractable 
in theoretical calculations. 

In principle, optical tomography by homodyne detec- 
tion can characterize the state of an optical field com- 
pletely [9j|. Working in this direction, Kim et al. 10] have 
proposed an experiment to obtain the fidelity directly by 
mixing the two fields to be distinguished at a beam split- 
ter and measuring the Wigner function of the output. 
The method is valid so long as one of the field states is 
pure. Our proposal, instead, is to perform independent 
homodyne measurements on the two fields, and deter- 
mine their closeness from the probability distributions of 
the measurement results. Homodyne detection is highly 
efficient compared with photon counting, and more im- 
portantly, it is genuinely a CV measurement, in the sense 
that it reveals the continuous (wave) nature of the field. 
We aim to determine to what extent such a measurement 
is optimal in distinguishing Gaussian states. In general, 
homodyne detection refers to any measurement schemes 
where the signal optical field is superposed with an aux- 
iliary mode that has the same central frequency as the 
signal field. In this paper, we consider the usual bal- 
anced homodyne detection (BHD) where a very strong 
local oscillator field is mixed with the signal. The quadra- 
ture amplitudes X^ = (ae~ l< ^ + a^e^)/2 are then mea- 
sured in this BHD, where a is the annihilation operator 
for the optical field and <f> is an adjustable phase. We 
compare the quantum fidelity with the minimum overlap 
Icf, = f dx^y/p^x^y/piix^), the minimum taken with 
respect to all measurement angles <f>. 

Twamley has calculated the fidelity of two undisplaced 
thermal states ^1|, an d Paraoanu and Scutaru that of 
displaced thermal states 0. Scutaru has obtained the 
fidelity of arbitrary single-mode Gaussian states . We 
consider the set of all pairs of single-mode Gaussian states 
that have identical values for the means of the quadra- 
ture amplitudes, i.e., (X < j > ) Pl = {X^) P2 . We identify the 
conditions under which homodyne detection makes the 
optimal distinction between these states. While the con- 
sidered set of states is restricted, we note that many 
CV quantum information protocols do not change the 
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mean amplitude of the field. Examples include (imper- 
fect) quantum teleportation |l4j] and Gaussian cloning 
each with unity gain. 
The paper is structured as follows. In Sec. II we briefly 
review the basic formalism of Gaussian states and the 
calculation of the fidelity for a general pair of Gaussian 
states. Homodyne detection is considered in Sec. Ill, 
where we construct the overlap of the probability distri- 
butions for Gaussian states of equal mean. The minimal 
overlap is compared with the quantum fidelity to iden- 
tify those cases in which homodyne detection is optimal 
in distinguishing the states. Pairs of radially symmetric 
(unsqueezed) Gaussian states with different means are 
also briefly considered. A summary and discussion is 
presented in Sec. IV. 

II. FIDELITY OF TWO SINGLE-MODE 
GAUSSIAN STATES 

A. Gaussian states 

A Gaussian state p is completely characterized by its 
first and the second moments. It possesses a Gaussian 
characteristic function, 

C(x) = tr{ P W(x)} = e *(R)* T e -k*T* T i (i) 

where R = (q, p) , with q and p position and momentum 
operators, respectively, a; is a real 2-dim row vector, and 
W(x) = e tR x is the Weyl operator. The canonical op- 
erators q and p are related to the quadrature amplitudes 
X = (a + a+)/2 = q/y/2 and Y = ~i{a - a + )/2 = p/y/2. 
The covariance matrix T is real, symmetric, and positive- 
definite, with matrix elements 

Tij = (ARiARj + ARjARi), = 1, 2. (2) 

The Heisenberg uncertainty relation must be incorpo- 
rated in order for T to represent a legitimate physi- 
cal state. Thus, with the canonical commutation re- 
lations represented by the symplectic matrix X, with 
[Ri,Rj] = iSy, i.e., 

the matrix T + zS (Ty -MSjj = 2(ARiARj}) is required 
to be positive semidefinite. The condition is necessary 
and sufficient for the covariance matrix T to represent a 
physical state [lq |. 

Consider now a real linear transformation S from the 
canonical operators R to another set of operators R! , 
with R! T = SR T . If the canonical commutation rela- 
tions are to be preserved ([R'^Rj] = iSy), S must sat- 
isfy the relation SYiS T = X. Such so-called symplectic 
transformations correspond to unitary Gaussian opera- 
tions generated by Hamiltonians quadratic in the opera- 
tors R 16]. The transformation S maps a real covariance 



matrix T to another such matrix, 

r-»r' = srs T . (4) 

If r + »S > 0, then r' + »E = S(T + iX)S T > 0; 
a physical covariance matrix remains physical under a 
symplectic transformation. 

Since five independent real parameters fully determine 
the characteristic function — three for the covariance ma- 
trix F and two for the mean values (J2) — it is not difficult 
to see that a general single-mode Gaussian state can be 
parametrically represented by a squeezed and displaced 
thermal state |17| . We may write the density operator in 
the form 

p = D(a)S(r,e)p T S\r,6)ft(a), (5) 

where 

n— v ' 

is a thermal state of mean photon number n, S (r, 9) = 
exp \[r/2){e %6 a' <2 — e~ l9 a 2 )] is the squeezing operator, 

and D{a) = exp (aa^ — a* a) is the displacement op- 
erator. The covariance matrix of the thermal state is 
diagonal — Ft = diag{7,7}, 7 = 2n + 1 — hence the co- 
variance matrix of the state p takes the factorized form 
[see Eq. gj] 

r = #s (7) 

where 

(cos0 -sin6»\ q-(V* \ (a\ 

^-^sin^ cos0 )' b -{0 (8) 

and the parameters 9 and s = e 2r denote the direction 
and degree of squeezing, respectively. The displacement 
operator leaves the covariance matrix unchanged, but 
fixes the mean field amplitude as (a) = a = a x + ia y . 
Thus, the general single-mode Gaussian state is param- 
eterized by the real numbers {7, s, 6, a x , a y }. The pu- 
rity of the state is determined solely by 7, with detT = 
7 2 > 1; a pure state corresponds to the case 7 = 1. 
Without loss of generality, the squeezing parameter can 
be restricted to s > 1. A state is squeezed whenever 
7/s < 1. 

B. Fidelity 

Given two quantum states, p\ and P2, the fidelity F is 
defined by 

F(pi,p 2 ) = te^y/plp2y/pl. (9) 

The fidelity F is continuous with respect to p\ and 
P2, and concave, i.e., F (Ej^/jj, a) > Y H p i F{p i ,a) with 
YuPi — 1. It has the following properties as a measure of 
the closeness of the two states: 



3 



(i) F(pi,p 2 ) = 1 if and only if p\ — p 2 ; more generally, 
< F{ Pl ,p 2 ) < 1. 

(ii) F(px,p 2 ) is symmetric, i.e., F(px,p 2 ) = F(p 2 ,pi). 



(iii) F(pi,p 2 ) = (v^ilpal^i) when one of the states is 
pure — i.e., when p\ = |^i)(^x|- 

(iv) F{pi, p 2 ) does not change under a unitary transfor- 
mation U—i.e., F(Up 1 W,Up 2 U i ) = F(p!,p 2 ). 

Although the fidelity itself is not a metric, the angle be- 
tween two states defined by A(pi,p 2 ) = cos _1 [F(/9i, p 2 )] 
satisfies the triangular inequality A(p, a) < A(p, r) + 
A(r,a). 

The main difficulty in calculating the fidelity comes 
from the square root of operators. For Gaussian states 
the difficulty is readily resolved, however, since the 
characteristic function of the square-root of a Gaus- 
sian state is also Gaussian. Specifically, by a successive 
use of the composition rule in the position representa- 
tion, (x\p!p 2 \y) = J dz(x\p 1 \z)(z\p 2 \y), with (x\p\y} = 
cxp[— (ax 2 + dy 2 + 2bxy) + Ix + ky + g], Scutaru 
showed that the quantum fidelity of two Gaussian states, 
with covariance matrices I\ and mean amplitudes <X; = 
(a ix ,a iy ) (i = 1, 2), is given by 



F 



where 



ex P [-/3 T (r 1 + r 2 )- 1 /3] , (10) 



A = det(ri + r 2 ), 
8 = (detTi - l)(detr 2 - 1), 
/3 = e*2 — ax. 



(11) 



It is also possible to obtain the fidelity by employing the 
Uhlmann theorem 3], which states that 



F(pi,p 2 ) = max { |v t ) i |0)}|(*| 



(12) 



where the maximum is taken over all possible purifica- 
tions, and \(j)), of states pi and p 2 , respectively, in 
an extended Hilbert space Hr ® H$ — i.e., where p\ — 
tr^X*! and p 2 = tr fl |</>)(0|. 



III. DISTINGUISHING GAUSSIAN STATES BY 
HOMODYNE DETECTION 

When homodyne measurement of the quadrature am- 
plitude X^ = X cos (f> + Y sin <p is performed on a Gaus- 
sian state pi (parameterized by {^i, Si,9i,oti Xl otiy }) the 
probability distribution pi(x ( j > ) is given by 



Pi(x<t>) = / dx 4 , +7T / 2 W i (x < f > ,x < p +n/2 ) , 



(13) 



where Wi is the Wigner function 
1 

72 



d 2 XCi(X)e 



,3A*-/3*A 



(14) 



with Ci the characteristic function; in this expression /3 = 
Px + if3y is a complex variable, related to x^ via 



cos sin <p\ ( j3 x 
- smcf) cos cf) J yft y 



(15) 



and A — X x + i\ y is also complex, with [equivalent to 
Eq. Q] 

_ \a'-\' ai -±\ri\ T 



d(X) = tr{ Pl D(X)} = e 



(16) 



where A denotes the row vector (X y , —X x ) and D(X) is the 
displacement operator. From Eqs. (|13[) - l|16[) . we obtain 
the probability distribution 



Pi(x<t>) = 



■nB, 



exp 



\2 

J+aT" 1 sin 2 (0-00], 



(17) 



with 

B t = -fi[si cos 2 (0 
oti,<t> = OL ix cos 4> + a iy sin 4>. (18) 

Working from Eq. 1)17(1 we find that the overlap 1$ of 
the probability distributions of two Gaussian states, p\ 
and p 2 , is given by 



I<t> = J dx^yjpi (x,/, )\Jp 2 (x<f,) 

■{BxB 2 )* exp( - 



Bi + B 2 



1 



Bx+B 2 



Pi > (19) 



where (3^ = a 2 ^ — ot\^. Restricting ourselves, then, to 
states of the same mean (/3^, = 0), we arrive at the result 



/ 



B l 



f(x) 



V2X 1 /* 

v/T+F 



(20) 



and the fidelity [Eq. JHI] is 



F 



1/2 



(21) 
1,2), 



where, in terms of the parameters {7^, Si,9i} (i 
A = 7i 2 + 7l + hi72^(si,s 2 ,0), 
5 = (7x 2 - 1) (7l - 1) ■ (22) 
with 9 = 9 2 — 9i and 

D(si,s 2 , 9) = si + s 2+ - si_s 2 - cos 29, (23) 

where 

s i± = s i ±St 1 (* = 1,2). (24) 

Note that the function f(x) [Eq. pif ] is concave. There- 
fore in any interval x e [a, 6], the minimal value of f(x) 
is realized at one of the endpoints, x = a or x = b. Our 
task is to compare the minimal value of 1^ with the fi- 
delity (|21|l to determine under what conditions they are 
equal. To illustrate how the minimum is achieved, we 
first consider some simple cases. 
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A. si = 1 and S2 > 1 

In this case the Wigner distribution of the state p\ is 
a radially symmetric Gaussian [Fig. 1 (a)]. It is easy 
to see, then, that the minimum 1^ is achieved for a 
homodyne measurement in either the "unsqueezed" or 
"squeezed" direction of the state p 2 — i.e., for <j) = 9 2 or 
<f> = 9 2 + 7r/2 (note that squeezing in the strict sense re- 
quires 72/S2 < 1). The minimum is therefore achieved 
with either B 2 /Bi = (72S 2 )/7i or B 2 /Bi = 72/(7132)- 

(a) 71 > 72: In this case the minimum occurs for 
<t> = 6>2 + 7r/2 and B 2 jB\ — 72/(7132). Then, by compar- 
ing f(B 2 /Bi) with F, in a tedious but straightforward 
manner, it is found that the minimum is equal to F when 
(i) 71 = 72 = 1 (in which case both states are pure) or (ii) 
s 2 = (71 - 7r 1 )/(72 - 72~ 1 )> witn 7i> 72 > 1 (both states 
are mixed). The equality does not hold and homodyne 
detection is not optimal when one state is pure and the 
other is mixed (71 = 1 < 72). If 71 = 72 > 1, then s 2 
must be unity according to (ii), and the two states are 
identical. 

(b) 71 < 72: In this case the minimum occurs for 
<j) = 9 2 and B 2 /B\ = (72S2V71. Checking for equal- 
ity with F, in a similar fashion, we obtain the conditions 
(i) 7l = 72 = 1 and (ii) s 2 = (72 - 72^)/ (7i - 7f 1 ) I with 
71j72 > 1- 

In summary, homodyne detection is found to be opti- 
mal for distinguishing two Gaussian states when 
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71 = 72 = 1, (i) 
S2+ = r i2 (71,72 > 1), (ii) 



(25) 



FIG. 1: Parameterization for two Gaussian states with the 
same mean. Contours of the Wigner function are plotted: 
(a) the special case where state pi is radially symmetric in 
phase space (si = 1), (b) the general case. 



where 



Tl2 = 



_ 72 - 72 



7i - 7i 



71 - 7i 

72 - 7 2 ~ 



(26) 



Note that when two pure states are compared [case (i)], 
the equality 1$ = F holds for both phase angles, <fi — 9 2 
and (j) = 2 + it/2. 



Equation (|27(l has a solution if and only if the inequality 
A\ + A\ — A\ > is satisfied. In our case, this condition 
may be expressed as 



2T 2 -T£>( Sl ,s2,60 + 2<O, 



where 



(29) 



B. 



si > 1 and S2 > 1 



In this more general case [Fig. 1 (b)] it is not straight- 
forward to see which phase angle gives the minimal value 
of I,/,. Nevertheless, we can obtain the conditions under 
which the equality min^i./^ = F holds in the following 
way. 

The question is whether there exists an angle <fi for 
which 1^ — F. Note first, from Eq. 1)20(1 . that the equal- 
ity is equivalents B 2 /B 1 = [F- 2 ± (F~ 4 - I) 1/2 ] 2 , an 
equation that can be arranged in the form 



by using 



A x sin 2<t> + A 2 cos 2(f> + A 3 = 0, 

B 2 _ 72 (s 2+ + s 2 _ cos 2(0 - 6 2 )) 
B x 71 (sx+ + sx- cos 2(0 - 6 1 !))' 



(27) 



(28) 



T = (7i/72)[F" 2 ±(F- 4 -l)^] 



1/212 



(30) 



Importantly, inequality l|29ll has only three independent 
parameters, 71, 72, and D(sx, s 2 ,0). Thus, the condition 
for min{0}/0 = F may be expressed in terms of 71 ,72, 

and D(si, s 2 ,9); in other words, we need not treat the 
parameters si, s 2 and 9 independently. We can deduce 
general conditions by considering a few simple cases, such 
as 9 = and 9 = it/2, instead of solving the complicated 
inequality J2HJ). 

Let us first consider the case 9 = 6>2 — 8\ = — i.e., 
the two states are squeezed in the same direction. Since 
the extremal values of B 2 jB\ are then 72S2/71S1 and 
72S1/71S2, this case is similar to that in the previous 
subsection. Only condition (ii) of Sec. IIIA is changed, 
with s 2 replaced by s 2 / s\ in Eq. I|25|) . to obtain 



S2S1 1 



sis 2 1 



ri 2 



(31) 
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FIG. 3: Normalized difference of the probability overlap 1 4, 
measured in homodyne detection and the fidelity F when one 
state is pure and the other is mixed (71 = 1, 72 = 4), for 
si — 2 and 6 — tt/3. The plot shows that the equality 1^ — F 
does not hold for any S2. 



must be satisfied 0. The case is illustrated by Fig. 4. It 
turns out that condition l(33|) is very restrictive and the 
parameters satisfying it occupy a volume of measure zero 
in the 5-dimensional parameter space {71, 72, si, s 2 , 6}. 



FIG. 2: Normalized difference of the probability overlap 
measured in homodyne detection and the fidelity F in the 
case of pure states (71 =72 = 1), for si = 2 and 6 — tt/3: 
(a) as a function of S2 and <j>, (b) as a function of <j> for S2 = 1.5 
(solid), S2 = 2 (dotted), and S2 = 4 (dashed). It is seen that 
there always exists an angle <f> for which 1$ = F. 



Secondly, with 9 = 7r/2, we similarly have 72/71S1S2 < 
B2/B1 < 72S1S2/71. Thus, condition (ii) in Sec. IIIA is 
changed to 

sis 2 + (s^)- 1 =r 12 . (32) 

Note that, according to the definition J22J), the term 
s 2 s^ 1 + sis^ 1 in Eq. l(3i"|) is equal to ~.D(si, s 2 , 0), 
while the term sis 2 + (S1S2) -1 m Eq. 132(1 is equal to 
|D(si, S2, 7r/2). Thus, we deduce the general condition 
ensuring min^j/^ = i* 1 in the following way. 

(i) 71 =72 = 1: When one of the states is pure the 
other must be pure also in order to obtain equality. Then, 
homodyne detection provides optimal distinguishability 
without regard to the parameters {s\, s 2 , &i, 62}- This 
case is demonstrated in Fig. 2. Thus, two pure Gaussian 
states are always optimally distinguished by homodyne 
detection. Also, when one state is pure and the other 
not (71 = 1 < 72 = 1) equality between min^j/^ and F 
never holds, as illustrated by Fig. 3: one pure and one 
mixed Gaussian state are never optimally distinguished 
by homodyne detection. 

(ii) 71 > 1 and 72 > 1: In the case of two mixed states, 
the condition 

£>( Sl ,s 2) 0)=2r 12 (33) 



C. Gaussian states of different mean 

Finally, we comment on the case where the means of 
the two states are not the same. In this situation, the 
measurement direction (f> yielding the minimal value of 1$ 
is a function of the five parameters {si, s 2 , #2 — 61,0.2 — 
ai}. If we confine ourselves to pairs of Gaussian states 
that are radially symmetric in phase space (si = S2 = 
1), then it is easily shown that homodyne detection is 
optimal if and only if 71 =72; the variances of the two 
states must be the same. Examples include all pairs of 
coherent states. 



IV. SUMMARY AND DISCUSSION 

We have investigated the conditions under which 
homodyne detection can optimally distinguish, in an 
information-theoretic sense, two Gaussian states of the 
same mean. We found that two pure Gaussian states 
are always optimally distinguished by homodyne detec- 
tion. On the other hand, if one of the states is mixed, 
homodyne detection is optimal only for pairs of states 
satisfying Eq. ((33(1 . a condition satisfied within a region 
of measure zero in the parameter space. 

In general, when one of the states to be compared is 
pure, say p\ = | v I / i)(\E'i|, an optimal POVM is provided 
by the two-component measurement Eq = |^ r i)(Wi|, 
Ei = I — I \l/ 1 ) 1 1 , for which it is easy to see that 

EL=o V / PiM\/P2(m) = V / (*ilM*i) = F. For exam- 
ple, when pi is the vacuum state, the measurement can be 
implemented with a photodetector, by discriminating be- 
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FIG. 4: Normalized difference of the probability overlap 1$ 
measured in homodyne detection and the fidelity F in the 
case of two mixed states (71 = 2, 72 =4), for si = 2 and 
9 — 7r/3: (a) as a function of S2 and (j>, (b) as a function of <f> 
for S2 = 1.4 (solid), S2 = 1.1 (dotted), and S2 = 3 (dashed). 
It is seen that there exists an angle <f> for which 1$ = F only 
when Eq. loot is satisfied, i.e., for S2 = 1.4. The inset expands 
the region around the minimum. 

tween no-click (Eq) and click (Ei) events. Such a POVM 
is clearly state-dependent, however, and requires prior 
knowledge of the states involved. It is not, therefore, of 
general practical use. (Reference p| also constructs an 
optimal POVM that depends on the states to be distin- 
guished.) Our approach, on the other hand, considers 
a predetermined class of measurements, namely, homo- 
dyne measurements, and we show that the considered 
measurements can always optimally distinguish two pure 
Gaussian states; the prior knowledge is only that the two 
states are pure. Of course, our approach is only possible 



because the optimal POVM is not, in general, unique. 
In spite of its merit, the practical applicability of our re- 
sult is still rather restricted. For example, setting aside 
technical imperfections in homodyne measurement, both 
imperfect quantum teleportation and quantum cloning 
output a mixed state, even when the input state is pure. 

We have shown that two pure Gaussian states of the 
same mean can be optimally distinguished by a pure CV 
measurement^I| but we considered a restricted class of 
homodyne detection only, and the question remains open 
as to whether or not it is possible to define a generalized 
CV measurement to distinguish between mixed states. 
One possible generalization is the POVM measurement 
where each component is given by E a = ^U^\a)(a\U 
and J d 2 aE a = I. Here, U is an arbitrary unitary oper- 
ator and I a) is a coherent state. This measurement can 
be implemented as follows. First, the signal field p s is 
subject to the unitary evolution U and then mixed at a 
50:50 beam splitter with an auxiliary mode in the vac- 
uum state |0). Balanced homodyne detections are per- 
formed on the two output fields from the beam splitter, 
one for the Xq quadrature and the other for the 
quadrature (or, more generally, one for and the other 
for X,f J+ K/2)- The joint probability then corresponds to 
tr{E a p s } = ±(a\U Ps W\a). 

Now, if the unitary operation U is restricted to the 
symplectic transformations, the above POVM measure- 
ments correspond to projections into the pure Gaussian 
states [2(]]. One can consider, for example, the squeezing 
operations U — S(r,9), for which, in the limit r — ► 00, 
the measurement becomes that of homodyne detection 
of the quadrature amplitude Xg, the measurement con- 
sidered throughout this paper. On the other hand, in 
the opposite limit, r — > 0, it becomes the so called het- 
erodyne detection where the Q-distribution of the state 
is measured @. In this generalization, it would be in- 
teresting to prove or disprove the conjecture that the 
measurement scheme with r — > 00 is optimal among all 
S(r, ^-parameterized POVMs. Such considerations ap- 
pear of value to applications beyond the present limited 
context of the distinguishability of states and they are 
thus left for future work. 

This work was supported by the NSF under Grant No. 
PHY-0099576 and by the Marsden Fund of the RSNZ. 
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